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Near-Source  Observations  of  Quarry  Explosions 

Lane  R.  Johnson  and  Michael  A.  Leonard 

Center  for  Computational  Seismology,  Lawrence  Berkeley  Laboratory, 
and  Department  of  Geology  and  Geophysics, 
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ABSTRACT 

An  experiment  was  performed  in  which  broadband  waveform 
data  were  recorded  by  a  seismic  network  surrounding  a  series  of  8 
chemical  explosions  in  a  quarry.  The  experiment  is  described  and 
all  of  the  waveform  data  are  presented.  A  velocity  model  was 
obtained  from  the  travel  time  readings  and  this  was  used  to  esti¬ 
mate  moment  rate  tensors  for  5  of  the  events  which  were  single 
explosions.  The  results  indicate  simple  source  mechanisms  for 
these  single  explosions,  particularly  the  deeper  events.  The 
efficiency  with  which  chemical  energy  is  converted  to  seismic 
waves  shows  a  strong  dependence  upon  the  depth  of  the  source 
and  this  is  caused  primarily  by  a  decrease  in  the  comer  frequency 
of  the  source  time  function  as  the  source  depth  decreases.  When 
scaled  by  the  size  of  the  explosion,  the  results  for  these  small 
chemical  explosions  are  in  reasonable  agreement  with  similar 
results  for  much  larger  nuclear  explosions. 


Introduction 

This  report  describes  some  of  the  results  which  have  emerged  from  a  spe¬ 
cial  experiment  in  which  a  series  of  chemical  explosions  detonated  in  a  quarry 
were  recorded  with  a  near-source  network  of  seismometers.  One  of  the  objec¬ 
tives  is  to  compare  the  results  of  this  experiment  with  results  from  similar 
experiments  involving  nuclear  explosions  and  to  determine  if  certain  fundamen¬ 
tal  aspects  of  elastic  wave  generation  by  explosions  can  be  more  easily  inves'.- 
gated  with  small  controlled  chemical  explosions  than  with  larger  nuclear  explo¬ 
sions. 

The  recording  experiments  at  the  quarry  were  very  similar  to  experiments 
performed  in  the  past  for  nuclear  explosions  at  the  Nevada  Test  Site,  with  all 
distances  and  depths  scaled  down  by  a  factor  of  about  10  for  the  quarry  experi¬ 
ments.  This  series  of  controlled  experiments  has  several  po'ential  advantages 
over  similar  experiments  that  have  been  performed  at  thr  Nevada  Test  Site. 
The  same  network  was  used  to  record  a  series  of  different  explosions,  so  the 
possibility  exists  that  the  separate  effects  of  source,  propagation,  and  recording 
site  can  be  isolated.  Three  separate  explosions  were  detonated  at  different 
depths  in  the  same  drilled  hole,  thus  permitting  a  direct  comparison  of  the 


effects  of  source  depth. 

The  main  purpose  of  the  present  report  is  to  describe  the  experiment  and 
the  data  which  it  produced.  Some  preliminary  estimates  of  the  moment  tensors 
for  some  of  the  source  everts  are  also  presented.  So  far  the  analysis  has  been 
mainly  concerned  with  the  question  of  how  the  source  mechanism  of  an  explo¬ 
sion  is  affected  by  its  depth.  This  particular  set  of  experiments  is  well  suited  to 
an  investigation  of  this  affect. 

Experiment 

In  the  fall  of  1988  a  series  of  chemical  explosions  were  detonated  in  the 
Kaiser  Permanente  Quarry  near  Menlo  Park,  California  (37.32  deg  N,  122.11 
deg  W).  The  explosions  were  arranged  by  Dr.  Willie  Lee  of  the  U.S.  Geologi¬ 
cal  Survey.  The  U.S.  Geological  Survey  deployed  five  different  types  of  instru¬ 
ments  over  a  broad  area  within  and  around  the  quarry  and  recorded  a  large 
number  of  data  channels  (Lee  and  Gibbs,  1989).  The  University  of  California’s 
part  of  the  experiment  was  to  record  broadband  waveform  data  at  a  network  of 
stations  close  to  the  explosions  at  distances  within  a  few  depths-of-burial. 

The  experiment  consisted  of  ten  different  explosive  events,  including  sin¬ 
gle  explosions  at  various  depths  in  drilled  holes,  regular  quarry  blasts  consisting 
of  simultaneous  explosions  in  clusters  of  shallow  drilled  holes,  and  one  event 
which  was  a  combined  quarry  blast  and  single  explosion  in  a  drilled  hole  with 
both  fired  at  about  the  same  time.  Table  1  contains  information  on  the  loca¬ 
tions  and  sizes  of  the  explosions.  The  accuracy  of  the  origin  times  is  better 
than  0.03  sec.  The  locations  have  a  relative  accuracy  of  about  1  meter.  The 
elevation  and  depth  both  refer  to  the  middle  of  the  explosive  material,  with 
elevation  being  referenced  to  sea  level  and  the  depth  being  referenced  to  the 
surface  at  the  drilled  hole. 

Two  different  types  of  events  were  involved  in  the  experiment,  a  single 
explosion  in  a  drilled  hole  and  a  quarry  blast.  For  the  purposes  of  this  study  a 
quarry  blast  is  defined  as  a  cluster  of  explosions  in  separate  shallow  drilled 
holes  detonated  at  approximately  the  same  time.  Note  that  events  KQ1,  KQ3, 
and  KQ5  were  single  explosions  which  were  all  detonated  in  the  same  drilled 
hole.  Also  note  that  event  KQ10  was  a  combination  of  a  quarry  blast  and  a 
single  explosion  which  were  slightly  separated  in  both  time  and  space. 

The  recording  network  which  collected  the  data  used  in  this  study  con¬ 
sisted  of  1 1  triggered  digital  event  recorders.  The  sensors  were  three- 
component  force-balance  accelerometers.  The  data  were  digitized  at  200  sam¬ 
ples  per  second  with  a  resolution  of  12  bits  per  sample  and  the  complete  sys¬ 
tem  response  was  fiat  to  acceleration  between  frequencies  of  0.2  and  50  Hz. 

The  same  network  was  used  to  record  all  events  in  the  experiment.  The 
locations  of  the  stations  in  the  network  are  listed  in  Table  2  and  plotted  in  Fig¬ 
ure  1.  Also  shown  in  Figure  1  are  the  locations  of  the  explosions  of  Table  1, 
with  the  exception  of  events  KQfi  and  KQ7.  Both  of  these  latter  events  were 
located  outside  the  boundary  of  the  recording  network,  so  data  from  these 
events  were  not  included  in  the  study.  The  stations  of  the  network  were 
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distributed  more  or  less  uniformly  over  an  area  of  slightly  less  than  1  square 
km,  and  for  all  of  the  events  there  was  at  least  one  station  within  250  meters  of 
the  epicenter  which  produced  usable  data. 

Because  the  experiment  was  conducted  in  a  large  quarry,  the  topographic 
relief  was  considerable.  Table  2  shows  that  the  maximum  elevation  change 
across  the  network  was  about  150  meters.  All  of  the  explosions  used  in  this 
study  and  stations  UCB1  through  UCB6  of  the  network  were  located  within  the 
quarry  and  at  approximately  the  same  elevation,  with  event  KQ9  and  station 
UCB6  being  slight  exceptions,  having  smaller  elevations  than  the  rest.  In  con¬ 
trast  to  this,  stations  UCB7  through  UCB11  of  the  network  were  located  out¬ 
side  the  quarry  on  its  rim  and  had  elevations  which  average  about  70  meters 
greater  than  the  stations  within  the  quarry. 

Data 

As  mentioned  above,  only  the  data  which  were  collected  from  the  8  events 
which  were  located  within  the  recording  network  are  included  in  this  report. 
Due  to  a  variety  of  malfunctions,  usable  data  were  not  obtained  from  all  of  the 
recording  stations  for  all  of  these  8  events,  with  the  number  of  stations  produc¬ 
ing  usable  data  ranging  between  5  and  11.  A  total  of  198  usable  seismograms 
were  obtained  from  the  8  events. 

The  waveform  data  which  were  recorded  from  the  8  events  are  shown  in 
Figures  2  through  25.  Relative  locations  of  the  stations  with  respect  to  these  8 
events  can  be  found  in  Tables  3  through  10.  In  Tables  3-10  the  range  column 
denotes  horizontal  distance  from  epicenter  to  receiver,  while  the  distance 
column  denotes  slant  distance  from  hypocenter  to  receiver.  In  Figures  2-25  the 
seismograms  have  been  arranged  in  order  of  increasing  range  and  scaled  by 
multiplying  each  seismogram  by  its  range.  Zero  time  corresponds  roughly  to 
the  detonation  of  the  explosion.  Because  of  differences  in  elevation  of  the 
recording  stations,  an  ordering  in  terms  of  range  is  not  always  an  ordering  in 
terms  of  distance,  and  thus  first  arrival  times  are  not  always  monotonic  in 
range. 

In  preparing  Figures  2-25  and  Tables  3-10  the  horizontal  components  of 
acceleration  were  rotated  to  form  a  radial  component  (positive  away  from  the 
source)  and  a  transverse  component  (positive  clockwise  about  the  source  as 
viewed  from  above).  The  vertical  component  is  positive  up.  Also  included  in 
the  tables  are  the  maximum  accelerations  measured  from  the  waveforms. 

Several  general  observations  apply  to  the  entire  data  set.  The  first  arriving 
signal  always  has  a  sense  of  first  motion  which  is  up  and  away  from  the 
source,  which  is  what  one  expects  from  an  explosive  source  below  the  receiver. 
Note  that  this  is  even  true  for  station  UCB6,  which  is  at  a  lower  elevation  than 
events  KQ2,  KQ4,  KQ5,  KQ8,  and  KQlOa,  indicating  that  in  its  path  from 
source  to  receiver  the  first  arriving  energy  has  followed  a  curved  path  and 
penetrated  to  a  lower  elevation  than  that  of  the  receiver.  This  is  to  be  expected 
if  the  velocity  increases  with  depth.  In  almost  all  cases  the  first  arrival  has 
larger  amplitudes  on  the  vertical  than  the  radial  component,  and  this  pattern  is 
most  pronounced  for  the  deeper  events  and  for  the  stations  at  higher  elevations. 
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On  the  transverse  components  the  first  discernible  motion  begins  at  the  same 
time  as  on  the  vertical  and  radial  components,  but  in  general  it  is  more  emer¬ 
gent  and  the  direction  of  first  motion  changes  from  station  to  station.  However, 
before  the  first  cycle  is  completed  the  amplitude  on  the  transverse  component 
has  usually  grown  to  a  value  comparable  to  that  on  the  radial  component. 

The  fact  that  event  KQlu  consists  of  two  separate  sources  is  clearly  evi¬ 
dent  in  Figures  23-25.  The  waves  from  the  shallow  quarry  blast  are  much 
smaller  in  amplitude  than  those  from  the  deeper  single  explosion  which  begin 
about  0.15  sec  later.  This  is  in  spite  of  the  fact  that  the  total  amount  of  explo¬ 
sive  in  the  quarry  blast  is  much  larger  than  that  in  the  single  explosion  (see 
Table  1). 

Some  of  the  advantages  of  using  the  same  network  to  record  a  series  of 
different  sources  are  readily  evident  in  the  waveform  data.  For  instance,  for  all 
of  the  events  station  UCB3  seems  to  have  anomalously  large  amplitudes  com¬ 
pared  to  other  stations  at  about  the  same  distance.  This  observation  applies  to 
all  three  components  and  suggests  the  existence  of  local  features  near  this  sta¬ 
tion  which  causes  an  amplification  of  the  seismic  waves.  In  principle,  station 
amplification  factors  such  as  this  could  be  extracted  from  the  data  set  itself 
using  a  least-squares  procedure.  Knowledge  of  such  factors  could  be  very  use¬ 
ful  in  the  analysis  of  the  data. 

The  data  set  is  particularly  well  suited  for  studies  of  how  the  depth  of  a 
source  affects  the  elastic  waves  it  radiates.  This  is  because  the  data  set  con¬ 
tains  records  from  3  different  explosions  of  about  the  same  size  detonated  at 
different  depths  in  the  same  drilled  hole:  KQ01,  KQ03,  and  KQ05.  The  depths 
(and  sizes)  of  the  explosions  were  217  meters  (1000  pounds),  106  meters  (900 
pounds),  and  42  meters  (900  pounds)  (see  Table  1).  A  comparison  of  the 
accelerations  from  the  deep  and  intermediate  explosions  (compare  Figure  2  with 
8,  Figure  3  with  9,  ant,  Figure  4  with  10)  shows  few  differences  with  similar 
waveforms  and  slightly  reduced  amplitudes  for  the  intermediate  depth  explo¬ 
sion.  However,  when  the  deep  or  intermediate  event  is  compared  with  the  shal¬ 
low  event  (compare  Figure  2  or  8  with  14,  Figure  3  or  9  with  15,  and  Figure  4 
or  10  with  16),  much  more  dramatic  differences  are  observed.  The  accelera¬ 
tions  from  the  shallow  event  are  much  reduced  in  amplitude  at  all  distances 
(note  the  factor  of  5  difference  in  the  scale  of  the  plots  for  the  shallow  event) 
and  the  seismograms  have  the  general  appearance  of  being  more  prolonged  in 
time. 

Velocity  Model 

The  seismograms  in  Figures  2-25  represent  the  combined  effects  of  both 
source  and  propagation.  One  of  the  objectives  of  this  study  is  isolate  the 
source  effect,  and  a  way  of  doing  this  is  to  solve  for  the  force-moment  tensor 
which  represents  the  source.  However,  this  requires  the  calculation  of  Green 
functions  and  this  in  turn  requires  that  a  velocity  model  be  available.  Thus  one 
of  the  first  steps  in  the  analysis  was  to  estimate  a  velocity  model  appropriate 
for  the  shallow  part  of  the  quarry  where  the  seismic  waves  were  propagating. 
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The  velocity  model  for  P-waves  was  estimated  from  the  travel  times  of  the 
first  arrivals  from  the  three  events  in  the  same  drilled  hole,  KQ01,  KQ03,  and 
KQ05.  These  travel  times  are  listed  in  Table  11.  The  inverse  problem  of  con¬ 
verting  these  travel  times  to  a  velocity  model  is  somewhat  different  from  that 
usually  encountered  in  seismology  in  that  the  large  variation  in  the  elevation  of 
the  receivers  must  be  taken  into  account.  The  problem  was  solved  by  using  a 
program  which  had  been  developed  for  estimating  velocity  structure  and  locat¬ 
ing  earthquakes  using  data  from  a  borehole  array.  The  problem  under  con¬ 
sideration  here  is  actually  easier  because  the  locations  of  the  sources  are  known 
and  can  be  fixed  at  the  outset.  The  basic  procedure  followed  by  the  program  is 
to  linearize  the  problem  and  perturb  the  velocity  through  a  series  of  damped 
iterations  until  the  travel  time  residuals  are  minimized  in  a  least  squares  sense. 
At  this  stage  of  the  analysis  it  was  decided  to  solve  for  the  average  one¬ 
dimensional  velocity  model  which  only  varies  in  the  vertical  direction. 

Figure  26  shows  the  one-dimensional  P-wave  velocity  model  for  the  lime¬ 
stone  rocks  of  the  quarry  in  the  vicinity  of  the  experiment  which  was  estimated 
using  procedure  described  above.  The  S  velocity  was  obtained  from  the  P 
velocity  by  assuming  a  constant  Poisson’s  ratio  of  0.25.  Note  that  the  shallow 
explosion  (denoted  by  the  number  5  in  Figure  26)  is  in  the  part  of  the  model 
where  there  are  steep  near-surface  velocity  gradients,  while  the  two  deeper 
explosions  are  in  a  part  of  the  model  where  the  velocity  gradients  are  more 
gentle.  In  this  respect,  the  shallow  event  is  similar  to  nuclear  explosions  at 
Pahute  Mesa  which  are  above  the  water  table,  while  the  two  deeper  events  are 
similar  to  events  which  are  below  the  water  table. 

Analysis 

So  far  the  data  from  this  experiment  have  been  analyzed  mainly  in  terms 
of  the  effect  which  source  depth  has  upon  the  explosion  mechanism.  The  basic 
idea  is  to  compare  sources  which  are  similar  except  for  their  depth.  While  a 
certain  amount  of  analysis  of  this  type  can  be  accomplished  by  comparing  the 
seismograms  themselves,  it  is  generally  more  useful  to  remove  the  effects  of 
propagation  in  order  to  allow  more  direct  comparisons  of  the  sources.  A  con¬ 
venient  way  of  performing  this  latter  process  is  to  estimate  the  force-moment 
tensor  of  the  source,  and  preliminary  results  of  this  type  have  been  calculated 
for  5  of  the  events. 

Given  seismograms  recorded  on  a  network  surrounding  the  source  and 
given  a  model  of  the  velocity  structure,  it  is  fairly  straightforward  to  estimate 
the  second-order  force  moment  tensor  of  the  source.  At  the  distances  involved 
in  this  experiment  the  seismograms  are  fairly  complex  and  contain  overlapping 
contributions  from  near-field  terms,  P-waves,  S-waves,  Rayleigh  waves,  and 
Love  waves.  Thus  the  Green  functions  used  in  the  moment  tensor  inversions 
must  be  fairly  complete  and  contain  all  of  these  different  parts  of  the  solution. 
These  can  be  obtained  for  one-dimensional  models  using  standard  propagator 
methods,  although  the  computations  are  quite  lengthly. 

The  considerable  variation  in  topography  across  the  recording  network 
presents  a  special  problem  for  the  calculation  of  the  Green  functions,  since  it 


-  6  - 


represents  a  departure  from  a  one-dimensional  structure.  A  compromise  solu¬ 
tion  to  this  problem  was  attempted  by  calculating  Green  functions  for  two 
different  models  of  the  quarry  velocity  structure,  one  having  a  free  surface  at 
an  elevation  of  360  meters  and  the  other  having  a  free  surface  at  an  elevation 
of  430  meters  (see  Figure  26).  The  first  model  was  used  for  the  stations  within 
the  quarry,  UCB1-6,  and  the  second  model  was  used  for  stations  outside  the 
quarry,  UCB7-11. 

In  representing  the  source  in  terms  of  its  second-order  force  moment  ten¬ 
sor  the  assumption  is  implied  that  the  dimensions  of  the  source  are  small  com¬ 
pared  to  the  wavelengths  of  interest.  For  a  frequency  of  50  Hz  the  S  waves 
used  in  this  study  have  wavelengths  of  about  20  meters.  For  the  events  in  the 
experiment  consisting  of  single  explosions,  the  assumption  that  the  source  is 
small  compared  to  this  length  seems  to  be  reasonable.  However,  for  the  events 
which  were  quarry  blasts  this  assumption  becomes  more  suspect.  It  is 
definitely  not  valid  for  event  KQ10  (see  Figure  1)  and  it  needs  further  study  for 
events  KQ02  and  KQ04.  For  this  reason,  this  report  contains  only  the  results 
of  moment  tensor  inversions  for  the  events  which  were  single  explosions. 

The  moment  rate  tensors  that  were  estimated  for  the  5  events  which  were 
single  explosions,  KQ01,  KQ03,  KQ05,  KQ08,  and  KQ09,  are  shown  in  Fig¬ 
ures  27-31,  respectively.  The  source  depths  for  these  events  ranged  between  10 
and  217  meters  and  the  explosion  size  ranged  between  300  and  1000  pounds. 

A  few  general  observations  apply  to  all  of  the  moment  tensor  results.  In 
all  cases  the  shape  and  polarity  of  the  signal  is  quite  similar  on  the  three  diago¬ 
nal  elements  of  the  moment  rate  tensor,  and  the  size  of  the  signal  on  these  ele¬ 
ments  is  larger  than  on  the  off-diagonal  elements.  This  result  implies  a  source 
mechanism  which  involves  a  spherically  symmetric  change  in  volume,  which 
agrees  with  the  simple  model  of  an  explosive  source  as  a  pressure  pulse  on  the 
interior  of  a  spherical  cavity.  The  two  deepest  sources  (KQ01  in  Figure  27  and 
KQ03  in  Figure  28)  have  particularly  simple  moment  rate  tensors,  consisting 
primarily  of  a  single  concentrated  one-sided  pulse  beginning  at  the  origin  time 
on  only  the  diagonal  elements.  This  corresponds  to  a  moment  tensor  consisting 
of  a  simple  step  having  a  rise  time  of  about  0.03  second.  The  moment  rate  ten¬ 
sor  for  the  shallower  events  are  more  complicated,  with  the  initial  pulse  still 
present  but  less  well  defined  on  the  diagonal  elements  at  the  origin  time,  and 
this  is  foilowed  by  other  longer  period  signals  which  are  generally  in  phase  on 
the  diagonal  elements  and  smaller  with  more  random  phase  on  the  off-diagonal 
elements. 

The  moment-rate  tensor  estimated  for  event  KQ08  (Figure  30)  is  the  most 
complicated.  The  initial  pulse  on  the  diagonal  elements  is  slightly  noncausal 
and  of  longer  duration  with  the  suggestion  that  it  may  consist  of  two  or  more 
pulses.  While  it  may  be  that  this  source  was  actually  more  complicated  than 
the  others,  there  are  other  possible  causes  that  must  be  considered.  For  the  pur¬ 
poses  of  the  moment  tensor  inversion  the  distribution  of  recording  stations  for 
this  event  was  much  worse  than  for  the  other  events,  with  only  15  components 
included  in  the  inversion  and  all  of  these  contributed  by  stations  located  in  the 
same  azimuthal  quadrant.  Thus,  this  moment-rate  tensor  is  judged  to  be  the 
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most  poorly  determined  of  the  5. 

A  simple  explosion  consisting  of  a  spherically  symmetric  change  in 
volume  should  have  a  moment  tensor  with  only  an  isotropic  part  and  a  zero 
deviatoric  part.  Here  the  isotropic  part  is  taken  to  be  one  third  the  trace  of  the 
moment  tensor.  As  pointed  out  above,  the  estimated  moment  rate  tensors  are 
all  dominated  by  their  diagonal  elements  and  the  signals  on  these  diagonal  ele¬ 
ments  are  quite  similar  in  both  their  shape  and  polarity.  This  means  that  the 
isotropic  part  of  the  moment  tensor  is  much  larger  than  the  deviatoric  part,  and 
this  is  consistent  with  expectations  for  a  simple  explosion.  Thus  the  isotropic 
parts  of  the  moment  rate  tensors  can  be  used  as  a  representation  of  the  explo¬ 
sive  part  of  the  source  mechanism 

The  isotropic  parts  of  the  moment  rate  tensors  shown  in  Figures  27-31 
were  calculated  and  their  Fourier  transforms  computed  to  obtain  the  spectra 
shown  in  Figures  32-36  for  the  events  KQ01,  KQ03,  KQ05,  KQ08,  and  KQ09, 
respectively.  All  of  these  spectra  are  similar  in  that  they  have  relatively  flat 
sections  at  low  frequencies  and  steep  decaying  sections  at  high  frequencies. 
The  low  frequency  levels  of  these  spectra  were  measured  in  the  5-10  Hz  range 
and  these  are  listed  in  Table  12.  The  high  frequency  asymptotes  were  extrapo¬ 
lated  back  to  where  they  intersected  the  low  frequency  levels  and  these  were 
used  as  estimates  of  the  comer  frequencies,  which  are  also  listed  in  Table  12. 
Also  shown  in  Table  12  are  the  maximum  values  of  the  isotropic  pans  of  the 
moment  rate  tensors  as  measured  in  the  time  domain,  with  this  maximum 
occurring  in  most  cases  during  the  pulse  which  begins  at  the  origin  time. 

The  spectra  of  isotropic  moment  rate  tensors  shown  in  Figures  32-36  and 
the  summary  values  in  Table  12  show  some  interesting  differences.  The  deeper 
events,  K()01  and  KQ03,  are  relatively  flat  out  to  comer  frequencies  between 
30  and  40  Hz.  In  fact,  the  comer  frequency  for  event  KQ01  is  close  enough  to 
the  comer  frequency  of  the  anti-alias  filter  at  50  Hz  so  that  part  of  the  roll  off 
at  high  frequencies  in  Figure  32  may  be  caused  by  the  recording  system.  The 
spectra  for  the  shallower  events  are  considerably  depleted  in  the  higher  frequen¬ 
cies  in  comparison  to  the  deeper  events.  For  event  KQ05  and  KQ08  the  comer 
frequencies  are  at  about  10  Hz.  Interpreting  the  spectra  of  event  KQ09  is  more 
problematical,  with  a  peak  in  the  spectra  near  10  Hz  and  what  appears  to  be  a 
comer  frequency  at  about  26  Hz.  It  should  be  emphasized  that  estimates  of 
spectral  characteristics  such  as  those  listed  in  Table  12  are  rather  subjective  and 
must  be  used  with  caution.  It  is  more  meaningful  to  make  direct  comparisons 
of  the  spectra  than  to  compare  estimates  of  the  spectral  characteristics. 

As  mentioned  above,  the  three  events  at  different  depths  in  the  same 
drilled  hole  (KQ01,  KQ03,  and  KQ05)  provide  a  well  controlled  test  of  the 
effects  of  source  depth  upon  an  explosive  source.  The  moment  rate  tensors  for 
these  three  events  (Figures  27,  28,  and  29)  exhibit  the  same  pattern  which  was 
observed  in  the  seismograms.  The  two  deeper  events  are  quite  similar  in  both 
the  shape  and  amplitude  of  the  time  dependence  of  the  moment  rate  tensors. 
However,  the  shallow  event  shows  significant  differences  from  the  other  two, 
with  smaller  maximum  amplitudes,  a  less  distinct  pulse  at  the  origin  time,  and 
relatively  more  energy  at  later  times.  Comparing  the  maximum  isotropic 
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moment  rate  tensors  in  the  time  domain  (Table  12)  shows  an  almost  linear 
dependence  upon  source  depth.  Recalling  that  the  size  of  the  three  explosions 
differ  by  only  10%,  these  results  indicate  that  the  efficiency  with  which  chemi¬ 
cal  energy  is  converted  to  seismic  waves  is  strongly  affected  by  source  depth. 

Additional  information  about  the  effects  of  source  depth  can  be  found  in 
the  spectra  of  the  isotropic  moment  rates.  Figure  37  compares  on  the  same  plot 
these  spectra  for  the  three  events  in  the  same  drilled  hole.  The  point  of  interest 
here  is  that,  while  the  three  explosions  have  similar  spectral  levels  for  frequen¬ 
cies  less  than  10  Hz,  they  differ  considerably  at  higher  frequencies.  In  particu¬ 
lar,  the  deep  and  intermediate  explosions  have  relatively  flat  spectra  out  to 
comer  frequencies  beyond  30  Hz,  while  the  spectrum  of  the  shallow  explosion 
is  markedly  deficient  n  higher  frequencies,  with  a  comer  frequency  which  is 
less  well  defined,  but  near  10  Hz.  However,  regardless  of  any  estimates  of 
corner  frequencies,  there  is  little  question  that  the  three  spectra  are  quite  similar 
below  10  Hz  but  show  marked  differences  above  10  Hz  which  appear  to  be 
related  to  source  depth.  These  results  are  entirely  consistent  with  the 
differences  that  were  noted  above  with  regard  to  maximum  accelerations  on  the 
seismograms  and  the  maximum  values  of  the  time  domain  isotropic  moment 
rates.  It  is  also  worth  noting  that  the  effects  of  source  depth  observed  in  this 
study  are  similar  to  those  observed  by  Flynn  and  Stump  (1988)  for  small  chem¬ 
ical  explosions  detonated  at  a  range  of  depths  in  dry  alluvium. 

In  this  study  the  method  of  moment  tensor  inversion  appears  to  have  been 
successful  in  removing  most  of  the  effects  of  wave  propagation  between  source 
and  receiver  arj  thus  allowing  a  more  direct  examination  of  the  source 
mechanism.  The  results  thus  obtained  seem  to  show  that  the  source  mechanism 
is  strongly  dependent  upon  source  depth.  However,  it  is  possible  that  the 
efficiency  of  the  inversion  pn>cess  could  depend  upon  the  depth  of  the  source 
in  such  a  way  that  differences  in  the  results  that  appear  to  be  related  to  the 
source  depth  may  have  been  introduced  in  the  process  of  doing  the  inversion. 
Such  reasoning  is  made  plausible  by  the  fact  that  the  accuracy  of  the  inversion 
process  is  directly  related  to  the  accuracy  of  the  Green  functions,  and  these 
may  be  less  accurate  for  shallow  sourees  where  scattering  from  near  surface 
three-dimensional  inhomogeneitics  is  likely  to  be  more  significant. 

An  independent  check  upon  the  validity  of  the  moment  tensor  results  can 
be  obtained  from  the  seismograms  themselves.  Figure  38  shows  the  first  arriv¬ 
ing  P  wave  at  the  same  station  for  the  three  explosions  in  the  same  drilled  hole. 
The  waveforms  have  similar  shapes,  although  the  pulse  from  the  shallow  event 
seems  to  be  composed  of  slightly  longer  periods.  Note  that  the  maximum 
amplitude  of  this  pulse  actually  decreases  as  the  distance  from  the  source 
decreases,  even  though  one  would  predict  the  opposite  effect  solely  on  the  basis 
of  geometrical  spreading  and  attenuation.  The  reasons  for  these  differences  are 
further  clarified  in  Figure  39,  which  compares  the  spectra  of  these  pulses  for 
the  shallow  and  deep  events.  The  spectrum  of  the  intermediate  event  is  very 
similar  to  that  of  the  deep  event.  Note  that  the  spectra  are  quite  similar  for  fre¬ 
quencies  below  10  Hz,  but  at  higher  frequencies  the  shallow  event  becomes 
increasingly  deficient  in  energy,  falling  off  by  a  factor  of  about  10  by  30  Hz. 
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These  results  are  almost  identical  to  those  shown  in  Figure  37,  and  thus 
confirm  the  general  pattern  observed  in  the  moment  rate  tensor. 

The  results  from  the  three  events  in  the  same  drilled  hole  suggest  a  couple 
of  simple  empirical  generalizations  about  the  effect  of  source  depth  upon  the 
explosive  source  mechanism.  First,  the  low  frequency  level  of  the  isotropic 
moment  rate  spectrum  is  unaffected  by  source  depth  and  is  related  to  the  size 
of  the  explosion.  Second,  the  comer  frequency  is  directly  related  to  the  source 
depth  and  decreases  as  the  source  depth  decreases. 

The  results  for  the  other  two  single  explosion  events,  KQ08  and  KQ09, 
seem  to  be  consistent  with  these  generalizations.  The  low  frequency  levels  of 
their  isotropic  moment  rate  tensors  are  similar  to  each  other  but  a  factor  of  2-3 
less  than  for  the  larger  events  KQ01,  KQ03,  and  KQ05  (see  Figures  35  and  36 
and  Table  12).  This  ratio  is  in  reasonable  agreement  with  the  relative  sizes  of 
the  explosions. 

The  comer  frequencies  of  the  isotropic  moment  rate  tensors  for  events 
KQ08  and  KQ09  are  also  consistent  with  the  results  for  the  other  three  single 
explosions  if  the  comparison  is  made  with  respect  to  the  relative  elevations  of 
the  events  rather  than  the  depths  below  the  surface  at  the  location  of  the  event. 
Thus  in  Table  12  event  KQ08  has  the  lowest  comer  frequency  and  the  largest 
elevation,  while  event  KQ09  has  an  elevation  in  between  KQ03  and  KQ05  and 
its  comer  frequency  also  falls  in  between  the  values  for  these  events.  This 
result  suggests  that  the  comer  frequency  is  more  affected  by  the  material  pro¬ 
perties  in  the  source  region  than  by  its  depth  below  the  surface.  If  this  were  so 
and  assuming  the  model  shown  in  Figure  26  is  reasonable  and  material  proper¬ 
ties  are  primarily  a  function  of  elevation,  then  one  would  expect  an  apparent 
dependence  upon  source  elevation  or  a  somewhat  more  ambiguous  apparent 
dependence  upon  source  depth. 

This  last  suggestion  that  material  properties  are  actually  the  controlling 
factor  in  determining  the  comer  frequencies  is  supported  by  the  observation  that 
events  KQ01,  KQ03,  and  KQ09  all  have  corner  frequencies  above  25  Hz  and 
are  all  located  at  lower  elevations  where  the  velocity  is  high  and  the  velocity 
gradient  is  low,  while  the  events  KQ05  and  KQ08  which  have  comer  frequen¬ 
cies  lower  than  10  Hz  are  located  at  the  higher  elevations  where  the  velocity  is 
low  and  the  velocity  gradient  is  high  (see  Figure  26).  It  should  be  noted  that 
for  sake  of  this  discussion,  ambient  stress  is  considered  to  be  a  material  pro¬ 
perty. 

Finally,  it  is  of  interest  to  compare  the  scaled  results  for  the  chemical 
explosions  of  this  study  with  similar  results  for  nuclear  explosions.  Such  a 
comparison  is  facilitated  by  converting  the  low  frequency  level  of  the  isotropic 
moment  rate  tensor  to  the  long  time  response  of  the  reduced  displacement 
potential  (see  Johnson,  1988,  for  details).  Then  letting  W  be  the  size  of  the 
explosion  in  kilotons  of  TNT,  the  results  shown  in  Figure  40  are  obtained. 
Note  that  in  order  to  have  a  consistent  depth  in  this  plot,  the  depths  used  for 
events  KQ08  and  KQ09  were  measured  from  the  top  of  the  drilled  hole  con¬ 
taining  the  other  three  events.  These  results  show  that  for  the  five  single  explo¬ 
sion  events  of  this  study  there  is  a  fairly  linear  relationship  between  the  scaled 
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long  time  response  of  the  reduced  displacement  potential  and  the  scaled  depth. 
According  to  simple  theory,  the  long  time  reduced  displacement  potential  is 
proportional  to  the  effective  volume  of  the  source.  Thus  these  results  indicate 
that  the  effective  source  volume  increases  as  the  source  depth  decreases.  This 
is  consistent  with  an  explanation  in  terms  of  decreased  overburden  pressure 
and/or  weaker  materials  at  shallower  depths. 

Also  shown  in  Figure  40  for  comparison  are  two  different  scaling  relation¬ 
ships  which  were  derived  from  nuclear  explosion  data.  The  scaling  relationship 
of  Mueller  and  Murphy  (1971)  is  partly  analytical  and  partly  empirical  and  is 
appropriate  for  large  nuclear  explosions  in  a  saturated  tuff-rhyolite  source 
medium.  The  Helmberger-Hadley  scaling  relationship  is  taken  from  Barker  et 
al.  (1985)  and  is  appropriate  for  the  same  type  of  source  medium.  An  interest¬ 
ing  aspect  of  Figure  40  is  that  the  scaled  results  from  this  study  exhibit  a  rela¬ 
tive  trend  that  agrees  well  with  both  of  the  scaling  relationships  developed  for 
nuclear  explosions  and  also  have  absolute  values  that  approximately  split  the 
difference  between  the  two  scaling  relationships. 

Figure  40  also  contains  scaled  results  from  two  nuclear  explosions  taken 
from  Johnson  (1988).  The  sizes  of  these  explosions  were  estimated  from  body 
wave  magnitudes  using  an  equation  from  Aki  et  al.  (1974).  This  gave  values 
of  107  Kt  for  Harzer  and  78  Kt  for  Chancellor.  The  values  were  taken 
from  Johnson  (1988)  and  were  obtained  by  the  same  type  of  moment  tensor 
inversions  used  in  this  study.  The  predictions  of  the  Mueller-Murphy  and 
Flelmberger-Hadley  scaling  relationships  are  also  shown  for  these  nuclear 
events.  Note  that  the  relative  agreement  between  the  scaled  moment  tensor 
results  and  the  published  scaling  relationships  is  about  the  same  for  both  the 
nuclear  explosions  and  the  chemical  explosions,  even  though  the  difference  in 
their  sizes  is  over  5  orders  of  magnitude.  It  seems  noteworthy  that,  although 
the  two  published  scaling  relationships  were  developed  specifically  to  explain 
large  nuclear  explosions,  they  do  almost  as  w’ell  in  explaining  the  small  chemi¬ 
cal  explosions.  This  seems  to  suggest  that  the  some  of  the  same  physical 
processes  are  involved  in  these  two  kinds  of  explosions.  Finally,  note  that,  in 
comparison  to  the  nuclear  explosions,  the  scaled  depths  of  the  single  explosions 
of  this  study  are  considerably  over  buried. 

Conclusions 

These  controlled  experiments  with  chemical  explosions  have  indicated 
that  source  depth  may  significantly  affect  the  efficiency  with  which  explosions 
generate  elastic  waves.  Furthermore,  this  effect  seems  to  be  strongly  frequency 
dependent.  The  results  generated  so  far  must  still  be  regarded  as  preliminary, 
but  they  have  suggested  several  interesting  questions  which  will  have  to  be 
examined  in  future  analysis.  A  critical  question  is  to  determine  what  physical 
processes  are  actually  controlling  this  apparent  source  depth  effect.  A  better 
understanding  of  these  processes  should  be  helpful  in  attributing  the  depth 
effect  to  the  relevant  physical  parameters,  such  as  overburden  pressure,  density, 
elasticity,  or  porosity. 
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An  impression  acquired  from  the  analysis  completed  so  far  is  that, 
with  respect  to  the  recorded  waveforms  and  the  estimated  moment  tensors, 
these  small  chemical  explosions  seem  to  exhibit  similar  characteristics  to  large 
nuclear  explosions.  However,  the  experiments  with  the  chemical  explosions  are 
much  easier  to  perform.  In  particular,  the  potential  for  isolating  and  studying  a 
particular  phenomenon  through  a  series  of  carefully  controlled  experiments  is 
much  higher  for  the  chemical  explosions.  Thus,  the  study  of  small  chemical 
explosions  offers  an  attractive  approach  to  investigating  some  of  the  remaining 
problems  associated  with  the  generation  of  elastic  waves  by  nuclear  explosions. 
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Table  1.  Shot  Information  for  the  1988  Kaiser  Quarry  Experiment 

Event 

Date 

OriginTime 

h:m:s 

Latitude 
deg  N 

Longitude 
deg  E 

Elevation 

m 

Depth 

m 

Explosive 

pounds 

KQ1 

Oct  27 

18:15:00.01 

37.32397 

122.10725 

136 

217 

1,000 

KQ2 

Oct  27 

18:30:16.42 

37.32182 

122.10890 

320 

10 

7,000/14  holes 

KQ3 

Oct  28 

18:15:00.01 

37.32397 

122.10725 

245 

106 

900 

KQ4 

Oct  28 

18:30:00.22 

37.32347 

122.10748 

335 

10 

1,900/4  holes 

KQ5 

Oct  31 

19:15:00.01 

37.32397 

122.10725 

311 

42 

900 

KQ6 

Oct  31 

19:30:00.16 

37.31998 

122.11247 

320 

10 

640 

KQ7 

Oct  31 

19:45:03.35 

37.32213 

122.12993 

525 

20 

1,510 

KQ8 

Oct  31 

19:46:00.16 

37.32377 

122.10752 

335 

10 

340 

KQ9 

Nov  1 

19:15:00.01 

37.320  5 

122.11037 

293 

12 

300 

KQlOa 

Nov  1 

19:30:00.81 

37.32408 

122.10752 

348 

10 

1 1,900/23  holes 

KQlOb 

19:30:00.91 

37.32380 

122.10707 

272 

81 

1,000 

Table  2.  Locations  of  Seismic  Recording  Stations 

Station 

Latitude 
deg  N 

Longitude 
deg  E 

Elevation 

m 

UCB1 

37.32357 

122.11045 

366 

UCB2 

37.32264 

122.11182 

366 

UCB3 

37.32190 

122.10810 

352 

UCB4 

37.32045 

122.10759 

357 

UCB5 

37.31860 

122.10723 

360 

UCB6 

37.32016 

122.11079 

306 

UCB7 

37.32351 

122.10526 

428 

UCB8 

37.32493 

122.10896 

4.35 

IJCB9 

37.32455 

122.11141 

454 

UC'BIO 

37.32115 

122.10524 

421 

UCB11 


37.32076 


122.10321 


404 
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Table  3. 

Maximum  Accelerations  from  Event  KQ1 

Station 

Azimuth 
deg  E  of  N 

Range 

m 

Distance 

m 

max  R 
cm/sec**2 

max  T 
cm/scc**2 

ir°*'  7 
cn  *2 

UCB7 

106 

183 

345 

134 

132 

329 

UCB8 

305 

186 

352 

115 

72 

264 

UCB3 

198 

241 

324 

197 

170 

393 

UCB1 

261 

287 

368 

53 

94 

112 

UCB10 

150 

359 

458 

35 

46 

140 

UCB9 

280 

374 

491 

59 

54 

185 

UCB4 

184 

392 

450 

37 

43 

117 

UCB2 

250 

431 

489 

70 

123 

137 

UCB11 

135 

505 

572 

26 

18 

62 

UCB6 

217 

526 

553 

25 

35 

34 

UCB5 

180 

596 

637 

31 

57 

79 

| 

Table  4. 

Maximum  Accelerations  from  Event  KQ2 

Station 

Azimuth 
deg  E  of  N 

Range 

m 

Distance 

m 

max  R 
cm/scc**2 

max  T 
cm/scc**2 

max  Z 
cm/scc**2 

UCB3 

82 

72 

79 

401 

287 

479 

UCB4 

143 

191 

195 

58 

58 

60 

UCB6 

222 

249 

249 

37 

41 

23 

L'CB2 

290 

275 

279 

17 

21 

39 

UCB8 

359 

345 

364 

19 

16 

11 

UCB9 

324 

376 

399 

7 

6 

8 

UCB5 

157 

387 

389 

16 

15 

23 

UCBII 

103 

518 

525 

9 

9 

9 

Table  5.  Maximum  Accelerations  from  Event  KQ3 


Station 

Azimuth 
deg  E  of  N 

Range 

m 

Distance 

m 

max  R 
cm/sec**2 

max  T 
cm/sec**2 

max  Z 
cm/sec**2 

UCB8 

305 

186 

266 

102 

66 

195 

UCB3 

198 

241 

264 

183 

200 

293 

UCB1 

261 

287 

311 

51 

63 

91 

UCB10 

150 

359 

400 

23 

32 

113 

UCB9 

280 

374 

428 

59 

41 

159 

UCB4 

184 

392 

408 

36 

42 

88 

UCB2 

250 

431 

448 

28 

31 

36 

UCBU 

135 

505 

529 

23 

14 

48 

UCB6 

217 

526 

529 

21 

39 

27 

UCB5 

180 

596 

607 

27 

42 

14 

Table  6.  Maximum  Accelerations  from  Event  KQ4 


Station 

Azimuth 
deg  E  of  N 

Range 

m 

Distance 

m 

max  R 
cm/scc**2 

max  T 
cm/sec*  *2 

max  Z 
cm/scc**2 

UCB3 

197 

182 

183 

94 

74 

106 

UCB8 

321 

208 

231 

21 

23 

22 

UCB1 

272 

263 

265 

17 

14 

10 

UCB10 

142 

324 

335 

14 

14 

26 

UCB4 

181 

335 

336 

20 

8 

21 

UCB11 

128 

483 

488 

8 

9 

7 
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l  able  7. 

Maximum  Accelerations  from  Event  KQ5 

Station 

Azimuth 
deg  E  of  N 

Range 

m 

Distance 

m 

max  R 
cm/scc**2 

max  T 
cm/scc**2 

max  Z 
cm/scc**2 

|  UCB8 

305 

186 

223 

55 

44 

72 

|  UCB3 

198 

241 

244 

64 

50 

90 

UCB10 

150 

359 

375 

13 

10 

31 

:  UCB9 

280 

374 

400 

9 

11 

14 

;  UCB4 

184 

392 

395 

16 

16 

27 

UCB2 

250 

431 

434 

8 

8 

11 

;  UCB11 

135 

505 

513 

6 

6 

15 

'  UCB6 

217 

526 

526 

8 

10 

10 

;  LCB5 

1 

1 _ 

180 

596 

598 

11 

18 

8 

Table  8.  Maximum  Accelerations  from  Event  KQ8 

Station 

Azimuth 
deg  E  of  N 

Range 

m 

Distance 

m 

max  R 
cm/scc**2 

max  T 
cm/scc**2 

max  Z 
cm/scc**2 

UCB3 

194 

213 

214 

35 

27 

43 

UCB10 

145 

353 

363 

10 

6 

8 

UCB4 

181 

367 

368 

8 

4 

6 

UCB6 

216 

494 

495 

3 

3 

2 

UCB11 

131 

507 

512 

4 

3 

3 
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Table  9.  Maximum  Accelerations  from  Event  KQ9 

Station 

Azimuth 
deg  E  of  N 

Range 

m 

Distance 

m 

max  R 
cm/scc**2 

max  T 
cm/sec*  *2 

max  Z 
cm/sec**2 

UCB4 

90 

246 

254 

22 

21 

42 

UCB3 

51 

258 

265 

72 

31 

102 

UCB2 

332 

275 

285 

72 

49 

72 

UCB5 

126 

346 

352 

27 

20 

12 

UCB10 

80 

461 

478 

8 

8 

21 

UCB9 

349 

464 

491 

14 

9 

17 

UCB11 

89 

635 

645 

6 

7 

19 

Event 

Elcv. 

m 

Depth 

m 

Explosive 

pounds 

Max.  Time  Domain 
10**20  dyne  cm/scc 

Low  Freq.  Level 
10**18  dyne  cm 

Comer  Freq. 
Hz 

KQ01 

136 

217 

1000 

1.56 

2.37 

37 

KQ03 

245 

106 

900 

0.81 

1.78 

32 

KQ05 

311 

42 

900 

0.31 

1.99 

11 

KQ08 

335 

10 

340 

0.09 

0.79 

10 

KQ09 

293 

12 

300 

0.28 

0.75 

26 
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Figure  1.  Map  showing  the  location  of  the  recording  stations  (UCB1  -  UCB1 1) 
and  the  location  of  the  explosions  (KQ1  -  KQ10)  at  the  Kaiser  Permanente 
Quarry. 
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Figure  2.  Radial  accelerations  of  ground  motion  for  event  KQ1  which  was  a 
1000  pound  explosion  at  a  depth  of  217  meters.  The  label  on  the  left  indicates 
the  recording  station.  The  accelerations  in  this  plot  have  been  scaled  by  multi¬ 
plying  by  the  epicentral  distance,  but  the  maximum  radial  accelerations  from 
this  event  were  0.20  g. 
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Figure  3.  Transverse  accelerations  of  ground  motion  for  event  KQ1  which 
which  was  a  1000  pound  explosion  at  a  depth  of  217  meters.  The  label  on  the 
left  indicates  the  recording  station.  The  accelerations  in  this  plot  have  been 
scaled  by  multiplying  by  the  epicentral  distance,  but  the  maximum  transverse 
accelerations  from  this  event  were  0.17  g. 


-21- 


UCB7 

l’CB8 

UCB3 

UCB 1 

UCB10 

UCB9 

UC.B4 

UCB2 


J _ I  V  I _ I  I _ L 


] 


50  cm/ sec  km 


0.0  0.1  0.2  0.3  0.4  0.5  0.6  0.7 

T i me ,  sec 


Figure  4.  Vertical  accelerations  of  ground  motion  for  event  KQ1  which  which 
was  a  1000  pound  explosion  at  a  depth  of  217  meters.  The  label  on  the  left 
indicates  the  recording  station.  The  accelerations  in  this  plot  have  been  scaled 
by  multiplying  by  the  epicentral  distance,  but  the  maximum  vertical  accelera¬ 
tions  from  this  event  were  0.40  g. 
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Figure  5.  Radial  accelerations  of  ground  motion  for  event  KQ2  which  was  a 
7000  pound  quarry  blast  distributed  in  14  holes  at  a  depth  of  10  meters.  The 
label  on  the  left  indicates  the  recording  station.  The  accelerations  in  this  plot 
have  been  scaled  by  multiplying  by  the  epicentral  distance,  but  the  maximum 
radial  accelerations  from  this  event  were  0.4 1  g 
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Figure  6.  Transverse  accelerations  of  ground  motion  for  event  KQ2  which  was 
a  7000  pound  quarry  blast  distributed  in  14  holes  at  a  depth  of  10  meters.  The 
label  on  the  left  indicates  the  recording  station.  The  accelerations  in  this  plot 
have  been  scaled  by  multiplying  by  the  epicentral  distance,  but  the  maximum 
transverse  accelerations  from  this  event  were  0.29  g. 
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Figure  7.  Vertical  accelerations  of  ground  motion  for  event  KQ2  which  was  a 
7000  pound  quarry  blast  distributed  in  14  holes  at  a  depth  of  10  meters.  The 
label  on  the  left  indicates  the  recording  station.  The  accelerations  in  this  plot 
have  been  scaled  by  multiplying  by  the  epicentral  distance,  but  the  maximum 
vertical  accelerations  from  this  event  were  0.49  g. 
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Figure  8.  Radial  accelerations  of  ground  motion  for  event  KQ3  which  was  a 
900  pound  explosion  at  a  depth  of  106  meters.  The  label  on  the  left  indicates 
the  recording  station.  The  accelerations  in  this  plot  have  been  scaled  by  multi¬ 
plying  by  the  epicentral  distance,  but  the  maximum  radial  accelerations  from 
this  event  were  0.19  g. 
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Figure  9.  Transverse  accelerations  of  ground  motion  for  event  KQ3  which  was 
a  900  pound  explosion  at  a  depth  of  106  meters.  The  label  on  the  left  indicates 
the  recording  station.  The  accelerations  in  this  plot  have  been  scaled  by  multi¬ 
plying  by  the  epicentral  distance,  but  the  maximum  transverse  accelerations 
from  this  event  were  0.20  g. 


Figure  10.  Vertical  accelerations  of  ground  motion  for  event  KQ3  which  was  a 
900  pound  explosion  at  a  depth  of  106  meters.  The  label  on  the  left  indicates 
the  recording  station.  The  accelerations  in  this  plot  have  been  scaled  by  multi¬ 
plying  by  the  epicentral  distance,  but  the  maximum  vertical  accelerations  from 
this  event  were  0.30  g. 
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Figure  11.  Radial  accelerations  of  ground  motion  for  event  KQ4  which  was  a 
1900  pound  quarry  blast  distributed  over  4  holes  at  a  depth  of  10  meters.  The 
label  on  the  left  indicates  the  recording  station.  The  accelerations  in  this  plot 
have  been  scaled  by  multiplying  by  the  epicentral  distance,  but  the  maximum 
radial  accelerations  from  this  event  were  0.10  g. 
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Figure  12.  Transverse  accelerations  of  ground  motion  for  event  KQ4  which 
was  a  1900  pound  quarry  blast  distributed  over  4  holes  at  a  depth  of  10  meters. 
The  label  on  the  left  indicates  the  recording  station.  The  accelerations  in  this 
plot  have  been  scaled  by  multiplying  by  the  epicentral  distance,  but  the  max¬ 
imum  transverse  accelerations  from  this  event  were  0.08  g. 
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Figure  13.  Vertical  accelerations  of  ground  motion  for  event  KQ4  which  was  a 
1 900  pound  quarry  blast  distributed  over  4  holes  at  a  depth  of  10  meters.  The 
label  on  the  left  indicates  the  recording  station.  The  accelerations  in  this  plot 
have  been  scaled  by  multiplying  by  the  epicentral  distance,  but  the  maximum 
vertical  accelerations  from  this  event  were  0.11  g. 
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Figure  14.  Radial  accelerations  of  ground  motion  for  event  KQ5  which  was  a 
900  pound  explosion  at  a  depth  of  42  meters.  The  label  on  the  left  indicates  the 
recording  station.  The  accelerations  in  this  plot  have  been  scaled  by  multiply¬ 
ing  by  the  epicentral  distance,  but  the  maximum  radial  accelerations  from  this 
event  were  0.07  g. 
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Figure  15.  Transverse  accelerations  of  ground  motion  for  event  KQ5  which 
was  a  900  pound  explosion  at  a  depth  of  42  meters.  The  label  on  the  left  indi¬ 
cates  the  recording  station.  The  accelerations  in  this  plot  have  been  scaled  by 
multiplying  by  the  epicentral  distance,  but  the  maximum  transverse  accelera¬ 
tions  from  this  event  were  0.05  g. 
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Figure  16.  Vertical  accelerations  of  ground  motion  for  event  KQ5  which  was  a 
900  pound  explosion  at  a  depth  of  42  meters.  The  label  on  the  left  indicates  the 
recording  station.  The  accelerations  in  this  plot  have  been  scaled  by  multiply¬ 
ing  by  the  epicentral  distance,  but  the  maximum  vertical  accelerations  from  this 
event  were  0.09  g. 
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Figure  17.  Radial  accelerations  of  ground  motion  for  event  KQ8  which  was  a 
340  pound  explosion  at  a  depth  of  10  meters.  The  label  on  the  left  indicates  the 
recording  station.  The  accelerations  in  this  plot  have  been  scaled  by  multiply¬ 
ing  by  the  epicentral  distance,  but  the  maximum  radial  accelerations  from  this 
event  were  0.04  g. 
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Figure  18.  Transverse  accelerations  of  ground  motion  for  event  KQ8  which 
was  a  340  pound  explosion  at  a  depth  of  10  meters.  The  label  on  the  left  indi¬ 
cates  the  recording  station.  The  accelerations  in  this  plot  have  been  scaled  by 
multiplying  by  the  epicentral  distance,  but  the  maximum  transverse  accelera 
tions  from  this  event  were  0.03  g. 
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Figure  19.  Vertical  accelerations  of  ground  motion  for  event  KQ8  which  was  a 
340  pound  explosion  at  a  depth  of  10  meters  The  label  on  the  left  indicates  the 
recording  station.  The  accelerations  in  this  plot  have  been  scaled  by  multiply¬ 
ing  by  the  epicentral  distance,  but  the  maximum  vertical  accelerations  from  this 
event  were  0.04  g. 
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Figure  20.  Radial  accelerations  of  ground  motion  for  event  KQ9  which  was  a 
300  pound  explosion  at  a  depth  of  12  meters.  The  label  on  the  left  indicates  the 
iccording  station.  The  accelerations  in  this  plot  have  been  scaled  by  multiply¬ 
ing  by  the  epicentral  distance,  but  the  maximum  radial  accelerations  from  this 
event  were  0.07  g. 
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Figure  21.  Transverse  accelerations  of  ground  motion  for  event  KQ9  which 
was  a  300  pound  explosion  at  a  depth  of  12  meters.  The  label  on  the  left  indi¬ 
cates  the  recording  station.  The  accelerations  in  this  plot  have  been  scaled  by 
multiplying  by  the  epicentral  distance,  but  the  maximum  transverse  accelera¬ 
tions  from  this  event  were  0.05  g. 
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Figure  22.  Vertical  accelerations  of  ground  motion  for  event  KQ9  which  was  a 
300  pound  explosion  at  a  depth  of  12  meters.  The  label  on  the  left  indicates  the 
recording  station.  The  accelerat’ons  in  this  plot  have  been  scaled  by  multiply¬ 
ing  by  the  epicentral  distance,  but  the  maximum  vertical  accelerations  from  this 
event  were  0. 10  g. 
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Figure  23.  Radial  accelerations  of  ground  motion  for  event  KQ10  which  was  a 
1 1,900  pound  quarry  blast  distributed  over  23  holes  at  a  depth  of  10  meters  fol¬ 
lowed  0.1  sec  later  by  a  1,000  pound  explosion  at  a  depth  of  81  meters.  The 
label  on  the  left  indicates  the  recording  station.  The  accelerations  in  this  plot 
have  been  scaled  by  multiplying  by  the  epicentral  distance,  but  the  maximum 
radial  accelerations  from  this  event  were  0.17  g. 
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Figure  24.  Transverse  accelerations  of  ground  motion  for  event  KQ10  which 
was  a  11,900  pound  quarry  blast  distributed  over  23  holes  at  a  depth  of  10 
meters  followed  0.1  sec  later  by  a  1,000  pound  explosion  at  a  depth  of  81 
meters.  The  label  on  the  left  indicates  the  recording  station.  The  accelerations 
in  this  plot  have  been  scaled  by  multiplying  by  the  epicentral  distance,  but  the 
maximum  transverse  accelerations  from  this  event  were  0.19  g. 
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Figure  25.  Vertical  accelerations  of  ground  motion  for  event  KQ10  which  was 
a  1 1,900  pound  quarry  blast  distributed  over  23  holes  at  a  depth  of  10  meters 
followed  0.1  sec  later  by  a  1,000  pound  explosion  at  a  depth  of  81  meters.  The 
label  on  the  left  indicates  the  recording  station.  The  accelerations  in  this  plot 
have  been  scaled  by  multiplying  by  the  epicentral  distance,  but  the  maximum 
vertical  accelerations  from  this  event  were  0.35  g. 
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Figure  37.  Comparison  of  the  amplitude  densities  of  the  isotropic  moment  rate 
tensors  for  the  three  events  in  the  same  drilled  hole,  KQ01,  KQ03,  and  KQ05. 
The  labels  on  the  curves  are  the  depths  of  the  events. 


Figure  38.  Comparison  of  the  first  arrival  on  the  vertical  component  at  station 
UCB3  for  the  three  events  KQ05,  KQ03,  and  KQ01.  The  depth  of  the  events 
are  shown  on  the  right.  The  pulses  have  been  scaled  to  have  the  same  ampli¬ 
tude  on  the  plot,  but  the  maximum  ground  accelerations  for  the  three  signals 
are,  from  the  top  down,  90,  293,  and  393  cm/sec**2. 


Figure  39.  Amplitude  density  spectra  of  the  pulses  in  Figure  38  for  the  events 
KQ01  and  KQ05. 
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Figure  40.  Scaled  long  time  level  of  the  reduced  displacement  potential  of  an 
explosion  versus  the  scaled  depth  of  the  explosion,  as  determined  from 
estimated  moment  tensors.  The  results  for  the  five  single  explosions  are 
labeled.  Shown  as  x’s  are  the  predictions  of  the  Mueller-Murphy  scaling  rela¬ 
tionship  and  as  o’s  are  the  predictions  of  the  Helmberger-Hadley  scaling  rela¬ 
tionship.  The  H  and  C  denote  results  for  the  nuclear  explosions  Harzer  and 
Chancel!  >r. 
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